BIORTHOGONAL ENSEMBLES 



Alexei Borodin 

Abstract. One object of interest in random matrix theory is a family of point 
ensembles (random point configurations) related to various systems of classical or- 
thogonal polynomials. The paper deals with a one-parametric deformation of these 
ensembles, which is defined in terms of the biorthogonal polynomials of Jacobi, La- 
guerre and Hermite type. 

Our main result is a series of explicit expressions for the correlation functions in 
the scaling limit (as the number of points goes to infinity). As in the classical case, 
the correlation functions have determinantal form. They are given by certain new 
kernels which are described in terms of the Wright's generalized Bessel function and 
can be viewed as a generalization of the well-known sine and Bessel kernels. 

In contrast to the conventional kernels, the new kernels are non— symmetric. How- 
ever, they possess other, rather surprising, symmetry properties. 

Our approach to finding the limit kernel also differs from the conventional one, 
because of lack of a simple explicit Christoffel— Darboux formula for the biorthogonal 
polynomials. 



1. Introduction 

Orthogonal polynomial ensembles are widely known. They are characterized by 
the property that the joint probability density of an ./V-point ensemble has the form 

N 

p(xi, . . . , x N ) = const ■ Y\_u(xi) Y\_(xi - x 3 ) 2 

i=l i<j 

for some positive weight function u(x) . The phase space / of an orthogonal ensemble 
is a finite or infinite interval of the real line. Usually the couple (w(x), I) corresponds 
to one of the classical systems of orthogonal polynomials. 

These ensembles play a very important role in the random matrix theory, see 
[Me]. They also serve as a rich source of various mathematical problems which in- 
clude Selberg integrals, differential equations for Fredholm determinants and many 
others, see, for example, [Me], [TW]. 

One of the main properties of orthogonal ensembles is the existence of simple 
formulas for all correlation functions. Namely, the nth correlation function has the 
form 

n 

(*) p nN (xi,... ,x n ) = Y\_^(xi) ■ dct[K N (xi,Xj)}" j=1 

2 = 1 
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where K^{x,y) is the Nth Christoffel-Darboux kernel; 

p N {x)p N -i{y) -PN-i{x)p N (y) 



K N {x,y) 



x-y 



{Pk{x)} is the system of orthogonal polynomials on I with weight tu(x) (see [D], 
[Me], [NW1]). Determinantal formulas (*) often allow to study the asymptotic 
behaviour of an TV-point biorthogonal ensemble when the number of points N goes 
to infinity. It turns out that after the appropriate scaling the correlation functions 
tend to a limit. The limit correlation functions also have determinantal form with 
a certain limit kernel. 

The well-known sine kernel 

sin7r(x — y) 
n(x - y) 

arises in the scaling limit of the classical polynomial ensembles in the bulk of spec- 
trum. In particular, the sine kernel arises in the case of the Hermite weight function, 
see, e.g., [Me, section 5.2]. 

At the left edge of the spectrum, both the Jacobi and Laguerre weights produce 
the same Bessel kernel ([F], [NW2] ) 

ipi{x)ip 2 (y) - <pi(y)<p 2 (x) 
x-y 

where 

<pi(x) = J a (2v/i) , ¥2 (x) = xip\ (x) ; 

a > — 1, J a {z) is the Bessel function. 

Both the sine and the Bessel kernels are symmetric and represent positive self- 
adjoint integral operators in 1? on R+ and on R respectively. 

In this paper we study a one parameter generalization of the orthogonal polyno- 
mial ensembles, the joint probability densities of our ensembles have the form 

N 

p h (x 1 , . . . ,x N ) = const Y[uj(xi)Y[ [(xi - Xj)(xf - x?)] 

i—l i<j 

where 9 is a fixed positive number. We call these ensembles biorthogonal. Orthog- 
onal polynomial ensembles correspond to 9 = 1. 

Biorthogonal ensembles inherit several nice properties from the orthogonal ones. 
For example, their correlation functions also have determinantal form, and corre- 
sponding kernels can be expressed via so-called biorthogonal polynomials, see [Kl] 
for general definitions. This fact was proved by K. A. Muttalib in [Mu]. He also 
argued that biorthogonal ensembles are of certain interest in physics. 

We consider three different cases: 

(1) J =(0,1), ^(x)=x Q , 

(2) /=(0,+oo), uj(x) =x a e- x , 

(3) J = (-00, +00), w(x) = \x\ a e~ x \ 
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Naturally, we shall call the second ensemble biorthogonal Laguerre ensemble, it 
depends on one real parameter a > — 1. 

Ensemble (1) is given by a special case of the Jacobi weight function. Namely, the 
factor (1 — x)P of the general Jacobi weight is absent in our formula for uj{x). Our 
techniques does not allow to handle the case of arbitrary 0, so we shall set = 
but we shall still use the words biorthogonal Jacobi ensemble in this case. Thus, 
our Jacobi biorthogonal ensemble also depends on one real parameter a > — 1. 

The weight function in (3) is, in contrary, more general than the Hermite weight 
function, which corresponds to a = 0. However, we shall deal with this more general 
case and use the words biorthogonal Hermite ensemble. This ensemble depends on 
one real parameter a > — 1. 

We prove that in all these cases after an appropriate scaling there exists a limit of 
the correlation functions, and we explicitly compute the limit kernels thus obtained. 

It turns out that the limit kernel in Jacobi and Laguerre cases is the same and 
it equals 



OO 

£ 

k,l=0 



(-lfx 



L,0l 



kW{ 



q+l+fc 



) lT(a + 1 + 61) a + 1 + k + 6V 



This kernel can be represented as the integral 

tf a '°\x,y) = 6 f i(st) • J a+1 . e ((yt) e )t a dt 

Jo 6 ' e 

(-x) m 



where 



J a ,b(x) = 



m=0 



m\T(a + bm) 



is the Wright's generalized Bessel function, sec [Wr], [E2, 18.1]. 
The limit correlation functions have the form 



pi a < a -V ( Xl ,...,x n ) = p L n a ^ ( Xl , . . . , x n ) = J] x? ■ det \lC^ ( Xi , x, ) 



i=i 



These functions are defined on R + . 

The limit kernel in the Hermite case is expressed via K^ a,e \x,y) and has the 
form 

The limit correlation functions are 

n 

p H n er{a > e) {x u . . . ,x n ) = TT |a;j| a • det \lC Her t a ' e \x i ,x j j 

They arc defined on the whole real axis. 

Two limit kernels \xy\ a / 2 IC Her( - a ^ (x,y) and {xy) a / 2 JC^ a - e \x, y) can be consid- 
ered as biorthogonal generalizations of sine and Bessel kernels, see below Examples 
5.5 and 3.5, respectively. 
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Though these kernels are not symmetric anymore, the transposition x <-> y leads 
to a non-trivial symmetry in the asymptotics of Laguerre and Hermitc ensembles. 
Namely, the change of parameters 

turns out to be equivalent to the transformation x i— > x e of the phase space, see 
Corollaries 4.7 and 5.6 below. Finite point ensembles do not possess this symmetry, 
it appears only in the asymptotics. 

The method that we use in this paper also provides a new approach to the 
asymptotics of classic orthogonal polynomial ensembles. 

The paper is organized as follows. Section 1 is the introduction. In Section 
2 we prove some general statements for further purposes. Section 3 deals with 
the biorthogonal Jacobi ensemble. In Section 4 we work with the biorthogonal 
Laguerre ensemble. The trick that we use there to obtain a formula for Christoffcl- 
Darboux kernels originates from our work on stochastic point processes arising in 
the representation theory of the infinite symmetric group (see [O], [B], [BO]). We 
hope to explain this connection in subsequent publications. In Section 5 we compute 
the asymptotics of the biorthogonal Hermite ensemble. Section 6 is an appendix, it 
contains rigorous proofs of two main theorems. Heuristic proofs of these theorems 
can be found at the appropriate places of the main text. 

I am very grateful to G. I. Olshanski and A. A. Kirillov for numerous and helpful 
discussions. 



2. Generalities 

We start with a discussion of well-known orthogonal polynomial ensembles, see 
[NW1], [Me]. Consider an ensemble of N points on a (possibly infinite) interval 
(a, b) of the real line with the joint probability density of the form 

N 

(2.1) p(xi, . . . ,x N ) = const ■ Y\_co(xi)Y[(xi - x 3 ) 2 

i—l i<j 



where w(x) is some positive weight function on (a, b). We can explicitly compute 
all correlation functions of this ensemble using orthogonal polynomials 

{pk(x)}kL , degpfe = fc; 

b 

J Pk(x)pi(x)uj(x)dx = 5 M - 

a 

Namely, we consider the Nth Christoffel-Darboux kernel 

N-l 

(2.2) K N (x,y)=J2Pi( x )Pi(y)- 

i=0 



The integral operator in L 2 ([a, b], uj(x)dx) with this kernel is the orthogonal 
projection on the A-dimensional subspace Span{l, x, . . . , x 1 ^^ 1 }. 
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b b 

(2.3) p n {xi, ■■■ ,x n )= ^fZTn\\ j ' " j ■ ■ ■ ,x N )dx n+l ■ --dx N 

a a 

have the following form (sec [D], [Me], [NW1]) 

n 

(2.4) p n {x u ... ,x n ) = J]_uj(xi) ■ det[K N {x i ,x j )] i y :j=1 . 

i=l 

The Christoffel-Darboux formula (see [S], [El]) 

Pn(x)pn-i(v) - Pn-i(x)pn(v) 



(2.5) K N (x,y) 



x-y 



provides a convenient way of analyzing the asymptotic behaviour of Kn{x, y) when 
N -> oo. 

Let us look at K^(x,y) from another point of view. Clearly, the kernel has the 
form 

N-l 

(2.6) K N {x,y)=Y J 

k,i=o 

where Cki are some constants. There is a simple way of saying what these constants 
are. Set 

b 



gtj = J x t+ ' J 'oj{x)dx. 



Let 

c = (°ki) k,rJo . G = (9ij ) ufJ ■ 

Proposition 2.1. With the preceding notation 
(2.7) C = G-\ 

Proof. By definition of the Christoffel-Darboux kernel, for all j — 0, 1, . . . , N — 1 

b 

J K N {x,y)y 3 dy = x 3 . 

a 

On the other hand 

6 



" " JV-1 N-l 

j K N (x,y)y : >dy = / ^ c kl x k y l ■ y J dy = ^ 



Ckigi 3 x k . 



k,l=0 k,l=0 
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By equating these two expressions, we get 

JV-l 

Ckigij = 5kj- □ 

1=0 

Proposition 2.1 and (2.6) provide another way of analyzing the kernel Kjy(x, y). 
As we shall see, this will also work for biorthogonal ensembles, where Christoffel- 
Darboux formula becomes rather complicated. 

Let us introduce the biorthogonal ensembles in the following very general way. 
Consider an ensemble of N points on (a, b) C M with the joint probability density 
of the form 

N 

(2.8) p b (x!, . . . ,x N ) = const ■ J| u(xi) ■ det[&(a:j)]^- =1 • det[r)i(xj)]^ j=1 , 

i=l 

where £,i(x), rji(x); i — 1,2, .. . are some functions defined on (a, b). The formula 
(2.1) is clearly a special case of (2.8), take £i(x) = r]i(x) = x 1 ^ 1 . 

Suppose, we managed to biorthogonalize and {77^} with respect to the pairing 

b 

(2.9) (tv) = J Z(x)ri{xMx)dx. 

a 

That is, we have two other systems of functions, say {Ci(^)}^i and {^i(x)}f =1 , 
such that 

Q E Span{£i, . . . ,&}, ipj e Span{i]!, . . . ,T]j}; 

= 8ij 

for all possible i and j. 

Exactly the same argument as for orthogonal ensembles proves that the correla- 
tion functions (defined as in (2.3)) of a biorthogonal ensemble have the form 

n 

(2.10) p b n ( Xl ,... ,x n ) = l[w(x i )-det[K b N (x i ,x j )}l j=1 , 

i=l 

where 

N 

(2.11) K b N {x,y)=Y,Ci(*)Mv)- 

is an analog of the Christoffel-Darboux kernel (2.2). 

In the next section we shall restrict ourselves to the case 

Zi(x)=x i - 1 , r H (x)=x ^; 6 > 0. 

Then {Q} and {if>i} are so-called biorthogonal polynomials, see [Kl]. In this case 
determinantal formula (2.10) was proved in [Mu]. 
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(In Section 5, however, it will be more convenient to take 

This transposition, clearly, does not change the ensemble (2.8) and its correla- 
tion functions (2.10), but it interchanges the variables x and y in the Christoffcl 
Darboux kernel (2.11).) 

Unfortunately, there is no simple analog of Christoffcl-Darboux formula in the 
biorthogonal case. For example, for integer 8 = k e N, the following formula is 
proved in [I] (cf. (2.5)) 

K n( x i V) = k Z k ( ^2 a rs(N+s+l{x)lpN-r{y) + (3( N (x)^ N+1 (y) 

x y \ 

\ r,s>0 
r+s<k— 1 

where a rs , (3 are some constants. But even in the simplest case 9 = 2 it requires a lot 
of technical work to compute these constants for known biorthogonal polynomials. 

We shall use another approach. The kernel K b N (x,y) has the following form (cf. 
(2.6)) 

N 

(2-12) K b N (x,y)= ]T c b Mx) m (y) 

k,l=l 

for some constants c b kl . Set 

b 

9% = J £,j{x)r)i{x)w{x)dx. 

a 

Let 

C = ( c ki)k,i=ii G = (g lj ) l:j = 1 - 
It turns out that we have the exact analog of Proposition 2.1. 
Proposition 2.2. With the preceding notation 

(2.13) C* = (G 6 )- 1 . 
The proof is just the same. 

To conclude this section, we shall consider even more general situation, where 
an analog of Propositions 2.1 and 2.2 holds. 

Let us fix two (possibly infinite) intervals {a\, b\) and (02, 62) of the real line. We 
consider a distribution p(xi, . . . , x^; 2/1, . . . , yjv) defined on (ai, bi) N x (02, 62)^ of 
the form 

(2.14) p(xi, . . . ,x N ;yi, . . . ,y N ) = const ■ det w(xi, yj) ■ det^(xj) ■ det r\ r i y y j ) , 

all subscripts vary from 1 to N. Here w(x, y) is a (generalized) function on (ai, 61) x 
(«2, ^2); Ci( a; )7 ^(a^) are some (generalized) functions defined on (a\,bi) and (02, b 2 ), 
respectively. 
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We choose the constant so that 

J p(x!, ... ,x N ;yi,... ,yjv)cfoi ■ --dx N dyi ■ --dy N = 1. 

(a 1 ,b 1 )«x(a 2 ,b 2 )« 

Suppose again that we managed to biorthogonalize the systems {&} and {rji} 
with respect to the pairing 

bi b 2 

(2.15) = j j ' €{x)r}(y)w(x,y)dxdy. 

ai a2 

Thus, we have somehow constructed two systems {Ci(x)}fL 1 and {ipi(y)}fL 1 , such 
that 

(i e Spanfa, ... , € Span{»7i, . . . , r]j}; 

(Ci.Vj) = <% 

for all possible i and j. 

If ai = a 2 = a, b\ = b 2 = b, and w(x,y) = w(x)S(x — y), then we return to the 
situation described above: (2.15) and (2.9) will coincide. 

It turns out that if we compute the "correlation functions" of the "measure" p 
only in (ai,6i) (i.e., we integrate p over all N y's and over some it's), then these 
correlation functions also have determinantal form. (In fact, as was recently proved 
in [EM], all correlation functions of measures of the type (2.14) have determinantal 
form, but this statement is much harder.) 

The following statement is proved in [MS]. 

Proposition 2.3. With the preceding notation 
N\ f 

^ N _ n ^ J p(xi, ■ ■ ■ ,x N ;yi, . . . ,y N )dx n+1 ■ ■ ■dx N dy 1 ■ ■ -dy N 

(a 1 ,b 1 ) N - n x(a.2,b 2 ) N 

(2.16) =det[K N (x hXj )}l j=1 , 
where 

N b2 

(2.17) K N {x,t)=J2Q(x) I \pi{y)w(t,y)dy. 
We shall use Proposition 2.3 in Section 4. 

Note that if a\ = a 2 = a, b\ = b 2 = b, and w(x, y) = uj(x)S(x — y), then 

K N (x,t) = Lu(t)K b N (x,t). 

As above, we can write 

N N 

(2-18) £Ci(aOVi(y) = Y, ^k(x)rn(y) 

i=l k,l=l 

for some constants Cki- 
Set 

bi 62 

9ij = J j ' £j(x)r)i(y)w(x,y)dxdy; 

c = (cki)k,i=i , G = (cjij )5=i ■ 
Quite similarly to Propositions 2.1 and 2.2 we get the following assertion. 
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Proposition 2.4. With the preceding notation 

(2.19) C = G-\ 
By (2.17) and (2.18) we get 

N b2 

(2.20) K N (x,t) = ~ c kitk{x) f m (y)w(t,y)dy. 

Note that everywhere above we supposed that the corresponding biorthogonal 
systems exist. The following simple statement will guarantee the existence in all 
our further examples. 

Proposition 2.5. // all principal minors of G are not zero, then there exist 
biorthogonal systems {Q{x)}f =1 and {^i{y)}iLi- 

Proof. The hypothesis implies that G possesses a Gauss decomposition: it can be 
represented as the product of a lower triangular and an upper triangular matrices. 
Thus, there exist a lower triangular matrix L = (lij) and an upper triangular matrix 
U = (uij) such that 

LGU = Id . 

Set 

N N 

^ = J2 u i^i ; '•'./ ^2 him ■ 
A straightforward check shows that these systems are biorthogonal. □ 

3. Biorthogonal Jacobi ensemble 

Our goal in this section is to study the asymptotic behaviour of the ./V-point 
ensemble on (0, 1) with the joint probability density (cf. (2.8)) 

N 



p J N ac (xi, ... ,x N ) = const ■ Y[ x? Y[ [{Xi ~ Xj){x e i - xf)] 

i—l i<j 
N 

= const j^J xf det x\~ x ■ det x ^ 

i=l 



where 8 > and a > — 1. We call this ensemble the biorthogonal Jacobi ensemble. 
Let us compute the matrix of pairwise scalar products (5^ ac ). We have 

l 

(3-2) ^ = 1*-™™*-**= . 1 . 

o 

To invert this matrix we shall use the following lemma. 
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Lemma 3.1. Let A = {A\, . . . , An} and B = {B\, . . . , Bn} be two sequences of 
complex numbers such that Ai + B 3 - ^ for all i,j = 1, . . . , N, and Ai ^ Aj, Bi ^ 
Bj for i 7^ j. Then 



1 ^ 1 



= (Cm) 



where 

(3.3) C u = 



A, + B } 

nr=i m+A^A+Bk) 



EM^ - Ai) W^kiBk - Bj) Ai + Bk' 

Proof. As is known, the elements of the inverse matrix are the ratios of the cofactors 
of corresponding elements of the initial matrix and the determinant of the initial 
matrix. The determinant of our matrix 

M 



is well-known Cauchy determinant, see [W]: 



det M 



UKA-A^-Bj)} 

N N ' 

n uiA+B,) 

z=ij=i 



Every submatrix of M has the same form as M itself for some other sets A and B. 
Then we can use the formula for the Cauchy determinant for computing any minor 
of M. In particular, we can compute all cofactors, and, thus, the inverse matrix. 
The result is exactly (3.3). □ 

Proposition 3.2. The inverse of the Gram matrix (#^ ac ) has the form 
(3.4) 



,,„,., , _ u{ (^)iv m-i)+a+i) N \ N 

l) + al 

k,l=i 



(gin = e 



(k-l)\(N-k)\ {l-l)\{N-l)\ k + 6{l-l) + a 



where (a) m = a(a + 1) ■ ■ ■ (a + m — 1) stands for the Pohgammer symbol. 
Proof. Direct application of Lemma 3.1 for 

Ai = 6(i-l), B, = i + a; i = l,...,N. □ 

Proposition 3.3. The correlation functions of the biorthogonal N -point Jacobi 
ensemble have the form 

n 

p J n a N c ( Xl , C„) ][,•'/. det W(x if Xj ))l i=1 

i=l 

where 
(3.5) 

K^tr ,a = e V (^)^v^ 1 (fl(*-i) + q + W (, - 1) (-i) fc+ ' 

N ^ (k-l)\(N-k)\ (l-l)\(N-l)\ k + 9(1-1) + a 
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Proof. The claim follows from Proposition 2.2 and Proposition 3.2. The existence of 
the Christoffel-Darboux kernel K^ ac (x, y) is guaranteed by Proposition 2.5 because 
all minors of the matrix (gff 10 ) are nonzero. □ 

Now we are in a position to compute the asymptotics of our ensemble as N — > oo. 
We shall employ the following entire function introduced by E. M. Wright, see 
[Wr], [E2, 18.1(27)] 

~ (-x) m 
J a ,b(x) = }2 m! r( a + 6m) 

m— 

(our notation differs from that used in [Wr], [E2]). It is closely related to Mittag- 
Lcffler type functions, see [E2, 18.1]. 
Note that 

(— l) m x m +f 
X* Ja+1,1 = V , — T = Ja(2Vx) 

m\T(a + m + 1 

m=0 x ' 

where J a (x) is the Bessel function. 

Theorem 3.4. For all n = 1, 2, . . . t/iere eiiste ffte Kmit 



x ' i=l 



where 



k ,4" n kW{s±^±)l\T(a + l + ei)a+l + k + ei 
(3.6) fe '' =0 e 



= [ Jb±i I (art) • J Q+ i.e((yt) e )t a rft. 
Jo 8 



Heuristic proof. If we substitute the asymptotic relations 

(fc-l)(8 + l)+o+l 

V 6 )N JV 



(jv-A;)! r(^a) 

(6{l-l) + a) N N (i-i)(e+i)+»+i 



{N-l)\ T{a + 1 + 6(1-1)) 

in (3.5) and shift the summation indices k and I by one, then we get (3.6). □ 

A rigorous proof of Theorem 3.4 (which is nothing more than a formalization of 
the argument above) can be found in Section 6. 

Example 3.5. For 0=1, Theorem 3.4 (and Theorem 4.5, see below) is well known, 
see [NW2] , [F] . In this case it is easy to check that 

x-y 

where 

Lpi(x) = J a (2y/x), lfi 2 (x) = XLp[(x), 
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is the Bessel kernel. 

Let us draw one more interesting corollary from Lemma 3.1. 

Note that in the notation of Section 2 if we know the Christoffel-Darboux kernels 
(2.11) then we can easily obtain explicit formulas for biorthogonal functions {(i(x)} 
and {ipi(x)}. Namely, we have the relation 

( m (x)il> m (y) = K h m {x, y) - K h m _ x [x, y). 

By Lemma 3.1 we can compute the kernel K^(x,y) for the systems 

£i(x)=x ai , r]i(x)=x bi ; 

a, + bj > -1, i,j = 1,2,... 

in L 2 ([0, 1], da;). The result of computation of the corresponding biorthogonal sys- 
tems is expressed in the next statement. 

Proposition 3.6. Let a\,a2, ■ ■ ■ and b\,b%, . . . be two sequences of complex num- 
bers such that 

a. + b^ > -1, i,j = 1,2, ... 
and ai ^ aj , bi ^ bj for i ^ j. Then two systems of functions 

n= 1,2,... 



ux) = Va n +b n + i± n^+^ + i) ^ 

and 

Mx) = Va n + b n + lJ2 I S? (bl t"" t \ ) x "^ " = 1 > 2 ,... 

f^i Uk=i,M bi ~ bk f 

are biorthonormal in L 2 ([0, 1], dx). In other words 

l 

J (m(x)lp n (x)dx = S mn . 


Example 3.7. For 

ai = bi = i - 1 

our biorthogonal systems are classic Jacobi polynomials. If 

a, = i — 1, bi = 0(i — 1) 



for some 9 > we get explicit formulas for biorthogonal Jacobi polynomials, cf. 
[MT1]. 
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4. BIORTHOGONAL LAGUERRE ENSEMBLE 

In this section we study the asymptotic behaviour of the ./V-point ensemble on 
(0, +oo) with the joint probability density 

N 

p L N a9 (x u ...,x N ) = const ■ TT [{ Xi - Xj ){xl - x e j)] TT< e-*i-- x » 

(4.1) ti ti 

= const ■ detxf 1 • detx" +e( ^ 1} • e~ Xl — - XN 

where a > — 1 and > 0. 

We call it the iV-point biorthogonal Laguerre ensemble. It turns out that the 
asymptotics of this ensemble is governed by the same kernel as that of biorthogonal 
Jacobi ensemble, see the previous section. However, in this case it is much more 
difficult to show. Here we do not know an analog of Lemma 3.1 for inverting 
the matrix of scalar products. However, the following trick allows to reduce the 
computation of the correlation functions to Lemma 3.1. 

Let us introduce N new variables j/i, . . . , t/jv an d a new distribution in 2N vari- 
ables 

N 

p N {xi,... ,x N ;yi,-.- ,vn) =p^ 9 (xi,... ,x N )Y[S(yi). 

i=i 

It is quite clear that 



<%n{xi, ... ,x n )= ^y^'^i J PN a9 ( x , ■ ■ ■ ,x N )dx n+1 ■ --dx N 



w 



N\ f 

_ n y J Pn{xi, ■ ■ ■ ,x N ;yi, . . . ,y N )dx n+1 ■ ■ -dx N dyx ■ --dy N . 



Thus, it suffices to compute the correlation functions of Pat. To do this we shall 
use Proposition 2.3, but, first, we need to show that pN has the form (2.14). 

Proposition 4.1. With the preceding notation 

p N {x\, . . . ,x N ;yi, . . . ,y N ) = const ■ det w(xi, yj) ■ det £i(xj) ■ det r)i(yj), 
all indices vary from 1 to N, and 

^{x)=x a+N+e ^- 1 \ r h {y) = 6«- 1 \v) 

for all i = 1, . . . , iV; 

g-X-J/ 

w{x,y) = — ■ — . 

x + y 

Proof. Using the formula for Cauchy determinant mentioned in the proof of Lemma 
3.1, we get 

II [{xi - Xj)(yi - yj)} n 

i<j - E {xi+yi) 

detw(Xi ' %) = — m^Ty-) — 6 

i,3 
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Next, 

N 

Y[(y i -y j ).det8^(y j ) = (-l)^N\l[6(y i ) 

i<j i—1 

because y m ■ 5^ n \y) = 8 mn for m> n. Finally, 



Pn(xi, ■ ■ ■ ,x N ;yi, . . . ,y N ) 



J ) ™ . , . N 



n(n-l) a+AT+0(*-l) *<j ' - E .TTj, . 

const- (-1) 2 detx v • e 1=1 • AT! I I <5(j/i) 

U( x i+yj) fJl 



hJ 

N N 



= const ■l[[(x i -x j )(x e i -x e j )] ■l[xf-e-^--^-l[5(y i ). □ 

i<j i—1 i—1 

Now, according to Section 2, we need to compute 

+■00 +00 

9ij = J f €j(x)rn(y)w(x,y)dxdy. 


Proposition 4.2. 

+ oc +00 

r r x a+N+6 ^-^5 { - i - 1 \y)e- x - y _ r(i + a + N + 6(j - 1)) 

' ' x~+~y V ~ a + N + 6(j-l)-(t-l) 





Before proving Proposition 4.2, let us explain our achievement. By Lemma 3.1 
we can now easily compute the inverse matrix 

C = (c fei ) = ((fa) -1 . 

Corollary 4.3. Mtt i/ie preceding notation 



(4.2) c w =0 



/ a+iV-a-l) \ 

rCl + a + ^fc-l))- 1 ( 5 ) N ^N+k+i+i 



(k - 1)\(N - k)\ {I - 1)!(JV — Z)! a + N + 6(k — 1) — (I — 1) 

Proof of Corollary 4. 3. Apply Lemma 3.1 for 

Ai = -{i-l), B l = a + N + i = l,...,N. a 

Proof of Proposition 4-2. It is easy to check that 

+oc +00 

f f x a v b e~ x ~ v 1 







r(a+l)F(6+l) x + y * a + 6+1 



if, say, 3?a and > 0. (Change of variables r = x + y, s = reduces the integral 
to the product of Euler's gamma and beta integrals.) 
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As is well known, there exists a distribution 



(j) c (u) 



r(c + i) 



which depends on c analytically, such that for c > — 1 it is just an integrable 
function u c /T(c + 1) for u > and for u < 0. For c < — 1 it is defined via analytic 
continuation. In particular, we always have the relation <p' c = 4> c -i and, thus, for 
any positive integer k 

Analytic continuation of (4.3) gives the relation 

+ OG +00 



X 



(y) — ; — dxdy = 







T(a+iy uyt " x + y " a + b+1 



for, say, 3?a > and 3fJ(a + 6 + 1) > 0. Our claim is the special case of this formula 
for 

a = a + N + 6(i-l), b=-j. □ 

Thus, we have inverted the matrix (g~ij), and as the result we get a formula for 
the correlation functions. 

Theorem 4.4. The correlation functions of the N-point biorthogonal Laguerre 
ensemble have the form 

n 

p L n a N 9 { Xl , ...,*„) = l[x?e-*< ■ Aet[K L N a9 { Xi , Xj )]l j=1 

i=l 

where 

Kjf(x, y) = E E T (a + 9k + l)r (*±§±I) 



fe!(iV — A; — l)!i!(r - i)! a + 6k + i + 1 ' 

Proof. By the definition of pjv, 

PnN( x l> ■■■ ,X n )= Tjy^-^i y ' PN ag ( x , ■ ■ ■ ,X N )dx n+1 ■ --dx N 

= 777 TT / P^v( x i> • ■ • ,x N ;yi,... ,y N )dx n+1 ■ ■■dx N dy 1 ■ --dy N . 

(N - n)\ J 

Thus, we can apply Propositions 2.3 and 2.4. We have 
y MvM*,v)dy = J S^Kyf-^-dy 



(-1) 



dy l 1 t + y 



,= hi 1 - 1 - 3 ) 1 



s-\„-t 
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Then, by Proposition 2.3, we get the determinantal formula 

PnN 9 ( X lT-- > X n) = det K N (Xi,Xj) 

where, see (4.2), 



N l-l 

k 



(/-!)! - • .* N 



k — X s — 



at i i , ( a+N-(l-l) \ 

A gr(l + a + 0(fc-l))- 1 ( 5 J 



JV 



,, , (fc-l)l(JV-fc)! (i-l-«)!(JV-0! 

( — 1)^ '+k+l + l x e(k-l)j.N -s-1 

X a + JV + 0(A:- 1)- (/- 1) ' 
Note that the factor (x/t) N disappears, when we take the determinant, and the 

n 

factor x a e~ t produces the factor J7 xf e~ Xi outside the determinant. Introducing 

i=l 

new summation indices 

i = N -I, r = N - s - 1 

and shifting index k by one, we arrive at our claim. □ 

One may ask whether the kernel K 1 ^ 9 (x,y) defined above is of the Christoffcl 
Darboux type in the sense of Section 2 (see (2.11)). The answer is positive: if we 
take (using the notation of Section 2) 

6(a) =!C *(*-i), r h (x)=x i - 1 

then 

N 

K L N ag {x,y) = YUx)Uy\ 



see Proposition 5.1. 

Now we are ready to pass to the limit N — > oo. 

Theorem 4.5. For all n = 1, 2, . . . i/iere exists the limit 



w/iere K^ a ' 6 \x,y) is defined in (3.6). 

Heuristic proof. Indeed, if we substitute in (4.4) the asymptotic relation 

{N-k-l)\ 
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then we get 

No \N<> J \N<> Ne J 

n—l , , >,-j.t Oh i n—1 



v ~ (-iy+ k x 9k y i y. 1 (y/my-* 

X k ^ k!r(a + 6k + l)i\T( l -±^)(a + 9k + i + l)f^ (r - »)! 

But 

and we arrive at our asymptotic kernel. □ 

A justification of the heuristic argument above is given in the appendix (Section 

6). 

Remark 4.6. Asymptotic behaviour of biorthogonal Jacobi and Laguerre ensem- 
bles is the same, while the scaling factors are different. 

Now let us point out a certain symmetry in the kernel K^ a,e \x, y) as defined by 
(3.6): 

It follows that the correlation functions (defined by the determinantal formula of 
Theorem 4.5), or rather the correlation measures 



■ dx\ . . . dx r , 



i=l 

are stable under the change of parameters 

(4.5) a ~^±i-i, e^- Q 

combined with the transformation of the phase space x i— ► x? , y i— » yi . 

For the Jacobi ensemble this is perfectly understandable, the transformation of 

i 

the phase space is equivalent to the change of parameters, because if X^ I ^ X J then 

N 

const Y\_ x ? TT _ X j)( X i ~ ' ' ' ^ X JV 

i— 1 i<_7 



const -JJ a;/ 1 JJ [(a* - ^)(a;/ - x? ) 
i=l i<j 



dxi • • • dxjy 



However, for the Laguerre ensemble this symmetry is rather surprising, the change 
of parameters is not equivalent to a transformation of the phase space, because the 
factor e -x\-—-x N dogg no t behave properly. Thus, we get a non-trivial conclusion, 
which deserves a separate statement. 

Corollary 4.7. The asymptotic behaviour of the Laguerre biorthogonal ensemble 
at the left edge of spectrum is invariant with respect to the change of parameters 
(4-5) and the transformation x i— » x~s of the phase space (0, +oo). 
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5. BlORTHOGONAL HERMITE ENSEMBLE 

Everywhere below the symbol x stands for sign(x)|a;| e , and for all integers k 

x ek = (x e ) k . 

In this section we are dealing with ensembles on (—00, +00) with the joint prob- 
ability densities 

N 



p% er {x u ...,x N )= const ■ JJ \ Xi \ a JJ [(*< - xj^xf - x*)] 

1 i<j 

N 

const ■ Y[ \xi\ a e~ Xi ■ detx^ 1 • dct 



(5.1) 

v N 



8(3-1) 
= 1 



where 9 > and a > — 1. We call these ensembles the biorthogonal Hermite 
ensembles. 

We shall reduce the study of their asymptotics to that of the Laguerre ensembles, 
see the previous section. To do this we need to introduce a notation for biorthogonal 
Laguerre and Hermite polynomials. 

Let us denote two sequences of biorthogonal Laguerre polynomials by {Z%(x, 9)} 
and {Y«(y,9)}. That is 

dcgZ:(x,9)=dcgY«(x,9)=n 

and 



p 00 

/ Z a m (x\ 0)Y°(x, 9)x a e- x dx = S mn . 
Jo 

Such polynomials were explicitly constructed in [K2], [C]: 



n\ (-lYx j 



w)=E0, rw+Q+1) . 



Note that our notation is slightly different from the conventional one, usually 
Z"(x,0) is multiplied by r ( en +^ +1 ) ; Du t we wan t our polynomials to be orthonor- 
mal. 

Following [MT2] we can construct biorthogonal Hermite polynomials {S%(x,9)} 
and {T°{y,9)} satisfying 



degS%(x,6)=degT%(y,9)=n 



/+00 
-OO 



S a m {x 6 ,9)T«{x,9)\x\ a e- x dx = 6 mn . 
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Namely, 

S$ n {x,e) = Z( a - 1 V\x 2 ,6); 
T? n (x,6) =Y^-V/\x\6); 
S^ n+1 (x,e)=xZ^/ 2 (x 2 ,e); 

Again, our notation differs from the usual one by scalar factors. 

It is quite clear that the Christoffcl-Darboux kernel for Hermite case can be 
expressed via that for Laguerre case, we shall do this in Proposition 5.3. 

Proposition 5.1. With the preceding notation 

N-l 

K^( x , y )=J2z?(x e ,9)Y«(y,6). 

i=0 

where K^ ag (x,y) is defined by (4-4)- 

Idea of the proof. One can prove this statement by direct verifying that 

K L N a ^{x,y)-K L N aa {x,y) 

is equal to 

we have explicit formulas for both these expressions. However, this check is rather 
tedious. □ 

Remark 5.2. Proposition 5.1 proves that K^ a9 (x,y) is the ChristofTel-Darboux 
kernel (2.11) for 

£.( x ) =x o(i-V, rh (x)=x i - 1 . 

Note that this fact and (2.10) provide another proof of Theorem 4.4. However, it is 
just a checking proof, the explicit formulas for biorthogonal Laguerre polynomials 
above do not prompt a suitable expression for the ChristofTel-Darboux kernel (4.4). 

In Section 4, using the trick with p N , we, actually, managed to derive the formula 
(4.4), and its relative simplicity allowed us to analyze the asymptotic behaviour of 
the biorthogonal Laguerre ensemble (Theorem 4.5). 

Moreover, we have a conceptual proof of the following statement which gener- 
alizes Proposition 5.1: let £i(x) = x ai for arbitrary complex numbers a^, 3?ai > 0, 
and r)i(x) = x' 1 ^ 1 ; then the trick with ppj described in Section 4 always produces 
the Christoffel-Darboux type kernel. But the proof is based on a certain formalism 
which exceeds the limits of the present paper. 

Set 

N-l 

K^ r (x,y)=J2sr(x ,0)Tt(y,e). 

i=0 

By (2.10), 

n 

p^( Xl , ...,x n )=H \ Xi \ a e- x ' ■ dct [K^ixuXi)]^ . 

i=l 
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Let us use a more detailed notation for Christoffel-Darboux kernels and write 

K L N a ^ a \x,y), K^ a \x,y) 

instead of 

indicating the dependence on a. 

Proposition 5.3. The Christoffel-Darboux kernel for the N-point biorthogonal 
Hermite ensemble has the form 

K^ a \x,y) = K L <^\x\y>)+x°yK L M aa ^\x\y*) for N = 2M 

< e, ' (a) {x, V) = K L <^\x\y^) + x e yK^ 9 ^\x 2 ,y 2 ) for N = 2M + 1 
where K^ 9 ^ a \x,y) is defined by (4-4)- 

Proof. Immediately follows from explicit expressions for biorthogonal Hermite poly- 
nomials via biorthogonal Laguerre polynomials. □ 

Now we can express the asymptotics of the Hermite ensemble via that of the 
Laguerre ensemble. 

Theorem 5.4. For all n = 1, 2, . . . there exists the limit (set M = ^ ) 

lim ^I^N*^^"^.'/. , 

N^oo M2S \M2S M20 J L J i,i=l 

where 

(5.2) K, Her ^ e \x,y) = K,(^< e ){x\y 2 ) + x e y^> e ){x\y 2 ) 

and K.( a > e \x, y) is defined in (3.6). 

The proof is straightforward. 

Example 5.5. The asymptotic kernel of the classic Hermite ensemble in the bulk 
of spectrum is the sine-kernel 

sin7r(£ — 77) 
tt(£ - 7]) 

Let us obtain it from our formulas. For = 1 we have, see (3.6), 

K^\x,y)= f J a+hl (xt)-J a+hl (yt)t a dt 
Jo 

But 

(5.3) JC Her ^(x, y) = JcH^V.y 2 ) + xy ■ JcH' 1 ) (x 2 , y 2 ) 
and 

00 ( x)^ 1 

- E fc!r(fc + 1/2) = cos{2 ^ 
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1 sin(2v^) 



k=0 



kT(k + 3/2) V% ' 



Thus, 



i 1 Z" 1 

/d - *' 1 )^ 2 ,^ 2 ) = - / cos(2xVt) cos(2yVi) 
t Jo 



eft 

1 ^sin2(.x — y) sin2(x + y) x 



2-7T \ x — y 



xyK.^ 1 ){x 2 ,y 2 ) = - I sm(2xVi)sm(2yVi)^= 

1 /sin2(a: — y) sin2(.T + y) 
2n \ x — y x + y 

Then (5.3) brings us to the sine-kernel for £ = 2x/ir, rj = 2y/ / K. 

Similarly to the Laguerre case, the biorthogonal Hcrmite ensemble also possesses 
a strange symmetry, cf. Corollary 4.7. 

Corollary 5.6. The asymptotic behaviour of the Hermite biorthogonal ensemble 
in the bulk of spectrum is invariant with respect to the change of parameters (4-5) 
and the transformation x i— ► sign(a;) • \x\~e of the phase space (—00, +oo). 

Proof. The claim easily follows from Corollary 4.7 and (5.2). If we set, see (4.5), 

a= — - 1 ' 0= 9 



then 

a— 1 i i i a-\-0 i i 



2 



These identities show that each summand of (5.2) is invariant under changes from 
the hypothesis. □ 

It would be very interesting to find some kind of natural explanation for Corol- 
laries 4.7 and 5.6. 

6. Appendix 

Proof of Theorem 3.4. The formula (3.5) for the Christoffcl-Darboux kernel 
K^ ac (x,y) implies that 

l 

(6.1) K J N ac {x,y) = 6 J A N (xt)B N {{yt) )t a dt 

o 

where 

N (k+a\ l_ T \k-l 
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(6-3) B N (y) = ^ { l- mN -iy. 



1=1 



Comparing (6.1) to (3.6) we see that it suffices to prove the following 
(6.4) lim — \+tA n ( * i ) = Jo+i i(x), 



^ j& n^ Bn {Jh) = Ja+lAy) 

where the convergence is uniform on every compact subset of M + . Indeed, these 
relations imply 

lim — ^ (*Xk^ (*,y\ = x «K<- a ' Xx,y), 

which is what we want to prove. 

We shall prove only (6.4); the proof of (6.5) is quite similar. 

Let us split the sum (6.2) into two parts: in the first part the summation index k 
runs from 1 to some M < N which will be chosen later, and the second part is the 
remainder. We shall denote these parts by A' N (x) and A'^(x), respectively. Thus, 

A N (x) = A' N (x)+A%(x). 

It will be sufficient to prove the following (uniform) estimates 

<"> -^M^)=l'l^£> +o(1,; 



< 6 - 7 > w-^G^tH 1 ' 

as N — > oo. 

To verify (6.6) we shall use Stirling formula 

(6.8) lnT(z) = (z-^Jlnz- m (27r) + o^ 1 ). 

Applying (6.8) to z — N + and toz = iV — fc + lwe get (we shall choose M so 
that it will be o(N), that is why the last term in (6.8) produces O(j^) in the next 
formula) 

r(iV+^) _( AT k + a l\,f AT k + a 



In i 2rf = \N + — In [N 

r(N-k + i)N h ^+ k - 1 V 92; 

-( N ~ k +^) ln (N - k + 1) - + k - lVln JV + 1) + O 
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Using asymptotic expansions 



ln(JV-fe+l) = ln7V r^ + O 



TV \ \ N 

we arrive at the following estimate 

T (N+^) k + a fk + a I 



In- 



T(N-k + l)N- 



i+fe-i N0 \ e 



n v 2y v ivy Viv/ Viv 



Now we want this expression to converge to as JV — > oo. Since k < M, we may 
set, for example, M = [N*]. Then 



and hence we get 



and 



N 



r (" + *t> = i + o ( ^) 

r(7v - fc + i)jv4 2 + fc - 1 



f M 1 V &k ( x \ k 

~ ^ r(7v - k + i)iv^+ fe - 1 (fc - i)!r 



M 



Thus, (6.6) is proved. 

To prove (6.7) we notice that for any a, b > 0, b is an integer < TV, we have the 
following simple estimate 

T(N + a) < (N + a) a+b+1 < ^ , a ^ «+*>+! 



r(7V-6)iV a + b - N a + b 



which can be obtained by applying applying the identity T(z + 1) = zT(z) to the 
numerator of the left-hand side [a] + b + 1 times if N — b > 1 , and [a] + b times if 
N-b = l. 
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Then 
1 



N 



- A" 



N 1+ v 

N 

* E 

fe=M+l 



AT 



v e In 



r(iv + ^) i^- 1 

r(N-k + i)7V^+fc-i (fc - i)ir (*£s) 



TV 

< e "(i 

fe=M+l 



AT0 



|fc-l 



(fc-i)ir(^) 



Since 1 + is bounded, say, by some constant c, and N < (M + l) 3 < fc 3 , the 
last sum does not exceed 



N 



E ^ 

fc=M+l 



- + fe 



lfc-1 



(fc-i)ir 

which is the difference of two partial sums Sn and Sm of the series 



E* 3 



ifc-i 



(fc-l)!r(4ja)- 

This series converges for all x uniformly on every compact set, and, consequently, 



N 



: a" 

x+i A N 



X 



N 1+ * 



converges to uniformly on every compact set, as was to be proved. □ 

Proof of Theorem 4.5. Let us rewrite the formula (4.4) for the kernel K^ ag {x, y) 
in the following form 

JV-l 



T(N) 



r (n + i±|+i; 



(6.6 



j + 1 2->. 



(r-i)! 



Note now, that if we substitute yTV 7 instead of y into the last sum, it will be close 
to 1: 

\ r — i 

V 



N-l 

y —J- 



(r - i)! ViV^ 



= 1 + 0(J\T*) 



where 0{N~~s) does not depend on i. Thus, we can neglect this sum while com- 
puting the limit. 

The rest of (6.8) can be written in the form (cf. (6.1)) 



JV-l 



*E 



T(N) 



r (at. 



i+a+l s 
, 



{-x e ) k {-yf 



(6.9) 



ki^o T ( a + dk + — k—l)\ T(N)T (i+s+i) i\ a + 9k + i + l 



= 9 / C N {(xt) e )D N {yt)t a dt 
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where 

N-l 



C N (x) = 



T(N)(-x) 



T(a + 6k + i)k\(N-k-l)V 

N—l 

D N (y) = £ 



N-l 



r 


[N + iH 


-a+l> 
8 / 


(" 




Y{N)T 


'i+a + 1 1 
v 8 y 


ii 



Comparing (6.9) to (3.6) we see that it suffices to prove the following 
(6.10) lim C N [A^ = J a+ i,i(x), 

JV— >oo \ JSje I 



(6.11) lim —^tDn [ JL ) = J«+i i (y) 

as N — > oo, because these relations imply the desired one 

Jim ^(^Vexpf^Wsf^-U =x^-%,*). 

N^oc No \NT> J \Ne J \Ne Ns J 

(The interchange x <-> y in the last expression does not change the determinants of 
the type det[K ( ~ a - 6 \x ll x j )].) 

The proofs of (6.10) and (6.11) are very similar to the proof of (6.4) which we 
carried out above, and we shall not give them here. □ 

References 

[B] A. Borodin, Point processes and the infinite symmetric group. Part IP. Higher correlation 
functions, Preprint, 1998. 

[BO] A. Borodin, G. Olshanski, Point processes and the infinite symmetric group. Part IIP. 
fermion point processes, Preprint, 1998. 

[C] L. Carlitz, A note on certain biorthogonal polynomials, Pacific Journal of Mathematics 
24 (1968), no. 3, 425-430. 

[D] F. J. Dyson, Statistical theory of the energy levels of complex systems I, II, III, J. Math. 
Phys. 3 (1962), 140-156, 157-165, 166-175. 

[El] A. Erdelyi (ed.), Higher transcendental functions, Vol. 1, Mc Graw-Hill, 1953. 
[E2] A. Erdelyi (ed.), Higher transcendental functions, Vol. 3, Mc Graw-Hill, 1953. 
[EB] B. Eynard and M. L. Mehta, Matrices coupled in a chain. I. Eigenvalue correlations, 



cond-mat/9710230 (1997). 
[F] P. J. Forrester, The spectrum edge of random matrix ensembles, Nuclear Physics B 402 

(1993), no. 3, 709-728. 

[I] M. N. Il'yasov, An analog of the Christoffel-Darboux formula for biorthogonal polynomi- 

als, Izvestiya Akademii Nauk Kazakhskoi SSR, Seriya Fiziko-Matematicheskaya (1983), 
no. 5, 61-64. (Russian) 

[Kl] J. D. E. Konhauser, Some properties of biorthogonal polynomials, Journal of Mathemat- 
ical Analysis and Applications 11 (1965), 242-260. 

[K2] J. D. E. Konhauser, Biorthogonal polynomials suggested by the Laguerre polynomials, 
Pacific Journal of Mathematics 21 (1967), no. 2, 303-314. 

[Me] M. L. Mehta, Random matrices, 2nd edition, Academic Press, New York, 1991. 

[MS] M. L. Mehta and P. Shukla, Two coupled matrices: eigenvalue correlations and spacing 
functions, Journal of Physics A - Mathematical and General 27 (1994), no. 23, 7793- 
7803. 



26 



ALEXEI BORODIN 



[MT1] H. C. Madhekar and N. K. Thakarc, Biorthogonal polynomials suggested by the Jacobi 
polynomials, Pacific Journal of Mathematics 100 (1982), no. 2, 417-424. 

[MT2] N. K. Thakare and H. C. Madhekar, Biorthogonal polynomials suggested by the Hermite 
polynomials, Indian Journal of Pure and Applied Mathematics 17 (1986), no. 8, 1031- 
1041. 

[Mu] K. A. Muttalib, Random matrix models with additional interactions, Journal of Physics 

A - Mathematical and General 28 (1995), no. 5, L159-L164. 
[NW1] T. Nagao, M. Wadati, Correlation functions of random matrix ensembles related to 

classical orthogonal polynomials, J. Phys. Soc. Japan 60 (1991), no. 10, 3298-3322. 
[NW2] T. Nagao, M. Wadati, Eigenvalue distribution of random matrices at the spectrum edge, 

J. Phys. Soc. Japan 62 (1993), no. 11, 3845-3856. 
[O] G. Olshanski, Point processes and the infinite symmetric group. Part I: the general 

formalism and the density function, Preprint, 1 998. 
[S] G. Szcgo, Orthogonal polynomials, 4th edition, American Mathematical Society, Collo- 

quium Publications, Volume XXIII, 1975. 
[TW] C. A. Tracy and H. Widom, Fredholm determinants, differential equations and matrix 

models, Comm. Math. Phys. 163 (1994), 33-72. 
[W] H. Weyl, The classical groups: their invariants and representations, Oxford University 

Press, 1939. 

[Wr] E. M. Wright, The asymptotic expansion of the generalized Bessel function, Proc. London 
Math. Soc. 38 (1934), 257-270. 

Department of Mathematics, The University of Pennsylvania, Philadelphia, PA 
19104-6395, U.S.A. E-mail address: borodine@matri.upenn.edu 



